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FOREWORD 


The analysis and development of the computer code described herein were 
initiated by W. R. Briley and H. McDonald while at the United Technologies 
Research Center, and completed by Briley and McDonald at the Scientific Research 
Associates Inc. under subcontract to United Technologies Research Center. In 
addition SRA Xnc. provided consulting assistance in the description of the com- 
puter code interfacing given in Section 9.7 and in its application to the test 
case described in Section 4.2. The details of the computing procedure with addi- 
tional results obtained by SRA are given in Section 10.0 together with the complete 
Sub-Contractor’s report. 
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1.0 SUMMARY 


An approximate analysis is presented for calculating three-dimensional, low 
Mach number, laminar viscous flows in curved passages with large secondary flows 
and corner boundary layers. The analysis is based on the decomposition of the 
overall velocity field into "inviscid" and "viscous" components with the overall 
velocity being determined from superposition. An incompressible vorticity transport 
equation is used to estimate inviscid secondary flow velocities to be used as 
corrections to the potential flow velocity field. A "parabolized" streamwise 
momentum equation coupled to an adiabatic energy equation and global continuity 
equation is used to obtain an approximate viscous correction to the pressure and 
longitudinal velocity fields. A collateral flow assumption is invoked to estimate 
the viscous correction to the transverse velocity fields . The approximate analysis 
is solved numerically using an implicit ADI solution for the "viscous" pressure 
and velocity fields and an iterative ADI procedure is used to solve for the inviscid 
secondary vorticity and velocity fields. 

This method was applied to computing the flow within (a) a turbine vane passage 
with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, 
and (b) a constant radius curved rectangular duct with R/D = 12 and 14 and with 
inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls. In the latter 
case a comparison is made with experimental data. 

A description of the operation and details of the computer code is also 
presented. 


1-1 



2.0 INTRODUCTION 


Future gas turbine engines will employ higher turbine inlet temperatures to 
achieve higher cycle efficiencies. An important consideration in the design of 
gas turbine engines therefore is the accurate prediction of loss and heat transfer 
in the region of the turbine end wall where large secondary flows and corner boundary 
layers occur. The flow in this region is three-dimensional and can not be treaded 
by an inviscid-viscous boundary layer interaction approach. Although direct 
numerical solution of the Navier-Stokes equations is possible, it is costly and 
time consuming on the computer. Therefore, efforts have been made to simplify 
the governing equations by neglecting streamwise diffusion and thus reduce the 
Navier-Stokes equations, which are elliptic, to a parabolic set of equations which 
can be integrated by forward marching numerical methods. 

Efforts to integrate these parabolic equations have been successful for 
straight ducts of rectangular cross section (Refs. 1, 2, 3). Even these methods, 
however, require further simplification of the problem by splitting the pressure 
into viscous and inviscid terms and treating the transverse momentum equations 
differently from the streamwise momentum equation. This method has been extended 
to ducts with other cross sections (Ref. 4) and to circular arc ducts of 
rectangular cross section with a large radius of curvature (Ref. 5). For all of 
these problems the secondary flow was small. 

In the turbine end wall problem the flow curvature and secondary flows are 
much larger and strong "elliptic" effects are present. Therefore it is to be 
anticipated that methods developed for flows with small curvature will encounter 
difficulty in these cases mainly due to the manner in which the local pressure 
distribution and transverse momentum laws are approximated. The purpose of this 
study is first to develop an approximate forward marching analysis for flows in 
ducts with large curvature and secondary flows and secondly, through numerical 
experimentation, to assess the analytic model and its capability for predicting 
end-wall losses and heat transfer in a vane passage. 
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3.0 ANALYSIS 


3.1 Basic Concepts 


A cross section of a simulated turbine vane passage duct geometry is shown 
in Fig. 1 and consists of two curved walls corresponding to the surfaces of the 
turbine airfoils. The remaining boundaries of the flow passage are two parallel 
flat surfaces comprising the endwalls. The geometry and inviscid flow, therefore, 
are two-dimensional. However, the viscous flow is three-dimensional, particularly 
near the endwalls, where the transverse pressure gradients associated with the ' 
turning of the flow produce strong secondary flow toward the low pressure or suction 
surface side of the passage. 


Coordinate System 


Curvilinear orthogonal coordinates x, y, z are constructed to fit the flow 

passage boundaries as shown in Fig, 1. Metric coefficients h]_, h 2 , I 13 are 

defined such that incremental distance 6 s is determined by (5s ) 2 = (hn 6 x ) 2 + 

0 0 

(l^dy ) 4 + (h^ 6 z) . In the planes parallel to the endwalls, orthogonal streamlines 
and velocity potential lines from a two-dimensional incompressible potential flow 
analysis are utilized as the coordinate lines for constant y and x, respectively. 
The z direction can be regarded as Cartesian, and thus hg = 1. The x coordinate 
is the primary flow or axial coordinate, and is associated with surfaces for 
which $■£ is constant. The y and z coordinates define transverse secondary flow 
planes at any given x location. In this coordinate system, the potential flow 
direction i^ and the normals i n to the transverse planes coincide with i, the 
unit vector in the x direction. A two-dimensional incompressible potential 
flow analysis and computer program developed by Anderson (Ref. 6 ) was employed 
in the present investigation, without modification, to compute the necessary 
coordinate data and potential flow solution. 


3.2 Governing Equations 

The complete derivation of the governing equation for the approximate flow 
field description employed in this study are presented in Section 10.0, Appendix A, 
of this report. The resulting equations as they apply to the present geometry are 
presented in detail in the following discussion. 

All variables used in the following governing equations are nondimensional 
having been normalized by the following reference quantities: distance, L r ; 

velocity, U r ; density, p r ; temperature, T r ; total enthalpy, U 2 ; pressure, p r U 2 
viscosity, p r . These quantities may be specified arbitrarily. This normaliza- 
tion leads to the following nondimensional parameters: Mach number, M; Reynolds 
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number) Re; Prandtl number, Pr; and specific heat ratio, y. These parameters 
are defined by 


M = U r /c , Re ^ r U r L r /^. r , Pr = c pMr /k , y = C p /C y (3>1) 

whe^e is the molecular viscosity, k is thermal conductivity, and c and c v 
are the specific heats at constant pressure and volume. The reference speed of 
sound, c, is defined by c^ = , where R is the gas constant, 

6 ^ o 

Since in the present application the coordinate lines coincide with the 
potential flow streamlines and potential surfaces, the following simplifications 
hold: Vj = Wj = u £ = 0. The velocity decomposition Eq. (10.3.2) can thus be expressed 
as 


0 = Tu +Jv + kw = 7( u x +u v ) + J(v s + V v ) + k(w s +w v ) 


(3.2) 


where i, j , k are unit vectors in the x, y, z coordinate directions, respectively. 
The incompressible streamwise vorticity equation, Eq. (10.2.6) can be written 
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This vorticity transport equation, Eq.' (3.3a), is derived in Section 10.0 with the 
following assumptions: 

1) the flow is incompressible; 

2) streamwise diffusion is negligible; 

3) the secondary vorticity vector lies in the direction of primary potential 
flow (i.e., coordinate direction through the same channel geometry; 

A) the principal streamline curvature is obtained from a primary potential flow 
solution; 

and is solved with the boundary condition £ = 0 on the wall. It should be noted 
that Eq. (3.3a) is an approximate representation of the full vorticity equation 
and as such is limited in its ability to properly model vorticity near the wall for 
two specific reasons. First, in the exact equation, Eq. (10. 2 . 6) , the convection 
velocities vanish at the wall, whereas in Eq. (3.3a) convection does not. Secondly, 
since vorticity can be generated only at a no slip surface, the boundary condition 
£ = 0 may not adequately model the generation of new vorticity at the wall. 
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In equation (3.3a) Q 2 


+ v 2 
I s 


■f w 
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and the vorticity is related to velocity by 
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where the scalar and vector potential functions <j> and <|» s are governed by 


and 
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(3.5) 

(3.6) 


Eq. 3.5 and 3.6 assume incompressible flow and are solved with slip boundary 
conditions. The streamwise momentum equation, Eq. (10.2.15) , is approximated by 
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Eq. (3.7) is solved with no slip boundary conditions and is derived in Section 10.0 
with the following assumption: 

1) Streamwise diffusion is negligible; 


2) the pressure can be split into viscous and inviscid terms in which the 
viscous term, p v , is a function of only the longitudinal coordinate, x. 
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An important result of this approximation is that transverse pressure gradients are 
set entirely by the calculated inviscid flow field pp and p g . 

The energy equation, Eq. (10.2.23) is approximated by 


Jr th ^ uE) + w [ | v ,(w * + V ' ()E ] + ir l h i h £^ (w s +w v ) E ] 
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Eq. (3.8) neglects streamwise diffusion of heat (heat conduction). It 
should be noted that a number of previous equations, namely Eq. (3.3) through 
Eq. (3.6), were derived for incompressible flow and therefore calculated results 
for temperature should be examined closely if compressibility effects are significant. 

The imposed pressure gradients are determined from an inviscid streamwise momentum 
equation, Eq. (10.2.18) given as 
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(3.9) 


Eq. .(3. .9.) is essentially the Euler equation for the primary potential and inviscid 
secondary flows and is related to the dependent variables using the following form 
of the auxiliary gas law relation, Eq. (10.2.25) 


a y-i u 2 \] ^Pj+P..) dp v (x) 

ax y P \ E ’ 2 /J = ax + dx 

The integral mass flux relation is 

Jf h g pudydz = constant 


(3.10) 


(3.11) 


To obtain the viscous contributions to the transverse velocity field, the 
momentum conservation law is dropped in favor of algebraic relations obtained 
from assuming that the flow is collateral (seeEq. (10.2.22), leading to the 
following equations for computing v v , w v ’ 


V 


v 



(3.12a) 
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3.3 Numerical Analysis 


The governing equations are replaced by finite-difference approximations. 
Three-point central difference formulas are used for all transverse spatial 
derivatives. An analytical coordinate transformation devised by Roberts (Ref. 7) 
is employed as a means of introducing a nonuniform grid in each transverse 
coordinate direction, to concentrate grid points in the wall boundary layer 
regions. Second-order accuracy for the transverse directions is rigorously main- 
tained. Two-point backward difference approximations are used for streamwise 
derivatives. The streamwise vorticity equation (3.3) is decoupled from other 
equations in the system and linearized with respect to £ by an ad hoc process 
consisting of lagging quantities not yet available at the implicit level. The 
resulting implicit difference equation for E, is solved using a scalar ADI scheme 
based on the technique of Douglas & Gunn (Ref. 8) for generating ADI schemes as 
perturbations of fundamental implicit schemes. Equations (3.5 - 3.6) for the 
scalar and vector potential functions <j> g and are elliptic in the transverse 
planes and are solved, given values for the right-hand sides, using scalar 
iterative ADI. Specifically, the Douglas-Gunn perturbation of the Crank-Nicolson 
scheme is used. The streamwise momentum (3.7), energy (3.8) and auxiliary 
gas law relation (3.10) are solved as a coupled system using linearized block ADI. 

A detailed discussion of the block ADI scheme used has been given by McDonald & 

Briley (Ref . 9 ) and Briley & McDonald (Ref . 10 ) . The general approach is to 

linearize the implicit equations by formal expansion about the solution at the 
most recent axial location. Terms in the difference equations are then grouped by 
coordinate direction, and the Douglas-Gunn (Ref. 8) technique is used to generate 
an ADI scheme. The resulting difference equations can be written in block-tridiagonal 
or a closely related matrix form and solved efficiently by block elimination techniques. 

A summary of the overall algorithm used to advance the solution a single axial 
step follows. It is assumed that the solution is known at the n level x n and is 
desired at x n+ ^-. 

1. The streamwise vorticity equation (3.3) is solved using scalar ADI to 
obtain 

2. The vector potential equation (3.5) is solved using scalar iterative ADI 
to obtain 
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3. yalue$ for v a and w s are computed using equation (3.4). Values of <j> s n 
obtained from the previous step are used. Values for v v and w v are computed 
from equations (3.12). The imposed pressure gradients 9(p-j- + p s )/3x are computed 
from equation (3.9). 

4. A value for the mean viscous pressure drop dp v (x)/dx is assumed. Initally, 
the value from the previous' step is used. 

5. Equations (3.7), (3.8), & (3.10) are solved as a coupled system using block 
ADI to obtain values of u n+ ^, p n+ ^ and E 11 ^. In general, the integral mass flux 
relation (3.11) will not be satisfied. 

6. Step 4 is repeated and this process is continued iteratively using the 
standard secant method (Ref. 11) to find the value of the mean viscous pressure drop 
which leads to u n+ -*- and p n+ '*' and satisfies the integral mass flux relation Eq. (3.11). 
The secant method was found in practice to converge to five figures on the third 
iteration. 

\ 

7. Using values now available for 9u/9x, the scalar potential equation (3.6) 

n4*1 

is solved using scalar iterative ADI to obtain <j> for use in the next axial 
step . 


To start this procedure, values for <f> n are needed. These are obtained by 
solving steps (4-7) above, under the assumption that the transverse velocity 
components are known and for these cases taken to be zero. 
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4.0 CALCULATION RESULTS 


The analysis for the flow field in a turbine vane passage in Section 3 must be 
considered approximate. A number of assumptions have been made including: (1) neglect 
of streamwise diffusion in the equations of motion, (2) decomposition of the velocity 
and pressure fields into viscous and inviscid terms, (3) use of an incompressible 
vorticity transport equation with inviscid velocities to estimate the secondary 
vorticity, (4) a' Viscous" static pressure field which varies only in 'the marching 
direction, (5) algebraic approximations to the’viscous" transverse velocity field. 
Several test cases are presented here to assist in assessing the range of applica- 
bility and utility of the approximate model used to determine the flow field 
properties . 

The first case is laminar flow through a simulated turbine vane passage; the 
second case is laminar flow through the same turbine vane passage with a higher 
inlet dynamic pressure. Finally, two cases are presented for laminar flow through 
a circular arc duct for which experimental data are available. 


4.1 Turbine Vane Passage 

The geometry of the simulated turbine vane passage is shown in Fig. 1. 

Upstream and downstream of the blades, the duct walls are assumed tangent to the 
inlet and exit velocities. The effect of the blade leading edge and trailing 
edge are neglected. Since the turbine geometric parameters such as twist, pitch, 
and chord are constant in the spanwise direction, the inviscid flow field is quasi 
two dimensional and can be constructed from the plane potential flow solution. 

This plain potential flow is shown on Fig. 1 and is used to construct an orthogonal 
coordinate system according to the method of Ref. 6 in which the stream function 
and velocity potential become the normal and streamwise coordinate. These stream- 
lines define the primary flow direction and the direction of the streamwise vorti- 
city. 

Small Inlet Dynamic Pressure 

A laminar flow calculation was performed on a 20 x 20 grid with 24 full 
streamwise steps or 48 ADI half steps. The initial conditions were for an inlet 
Mach number, M, = .1 and a Reynolds number. Re, = 1000. The remaining parameters 
for this test case are shown on Table I. As can be seen from Table I, the 
viscosity was adjusted to fix Reynolds numbers and Mach numbers and hence does 
not represent a value for air at these conditions. The run time for this case on a 
UNIVAC 1110 was 40 minutes. This run time was found to vary somewhat from case to 
case depending on the rate of convergence of the internal ADI iteration process. 

The calculated streamwise vorticity distribution at the last two stations is 
shown on Fig. 2. This vorticity is the component aligned in the primary flow direc- 
tion and produces the secondary flow velocities which lie in planes normal to the 
primary flow as defined in Refs. 12, 13. A large region of vorticity is apparent 
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on the suction side of the turbine endwall passage which may be identified as vortex 
roll up. The computational grid used in obtaining this result is superimposed on 
the solution at the last station. Grid points close to the wall are not shown 
because they are too dense to plot. These results appear qualitively reasonable 
within present knowledge of the flow field. 

The quantitative distribution of secondary vorticity is shown in Fig. 3. 

Again the abscissa is the distance across the gap and several slices are shown at 
increasing distance from the end wall out to the plane of symmetry. Note that 
the axis has been shifted vertically to assist in interpretation of results. In 
general it is seen that regions of vorticity occur near the suction side of the duct 
passage consistant with the vortex roll up process observed in Fig. 2. However, 
it is observed that the vorticity distribution is very irregular and in some cases 
two extrema appear with an intermediate reversal in sign for some cases. It would 
appear that the rather rapid decreases in vorticity, |, near the walls and subse- 
quent rapid increase across the gap is a direct result of the use of the inviscid 
boundary condition of zero vorticity on the walls - thereby raising some concern 
about the quantitative reliability of the approximate model being used to determine 
the local flow field properties. 

The gapwise and spanwise secondary flow velocities obtained using this vorti- 
city field are shown on Figs. 4 and 5, respectively . The secondary flow is seen to 
move across the gap from the pressure to suction side, then spanwise along the 
suction surface, back across the gap near the mid plane and finally down along the 
pressure surface. The grid points are shown as data points indicating sufficient 
resolution. Such a secondary flow field does not appear to be consistant with the 
secondary flow analysis developed in Refs. 12 and 13. 

The streamwise velocity distribution of the primary flow for this case is shown 
on Fig. 6. The abscissa is the distance across the gap, and several slices are 
shown at increasing distance from the end wall out to the plane of symmetry. 

These velocity distributions indicate that relatively thin boundary layers exist 
throughout the flow field despite the rather low Reynolds number of 1000. Also note 
should be taken of a small local extrema in the velocity distribution in the suction 
corner which may be caused by the large secondary vorticity. However, it should be 
noted that the viscous transverse velocities are not obtained by integrating the 
vorticity equation Eq. (3.3b) but rather by using an algebraic relation Eq.(3.12). 

In effect, the two transverse momentum equations have been replaced by two 
approximate algebraic relations. Since the transverse momentum equation relates 
pressure gradient, acceleration, and diffusion of vorticity; care should be taken 
in constructing this approximate model to show that the transverse momentum 
equations are satisfied. 

In an attempt to assess the internal consistency of the approximate model 
being used here, attention was given to its ability to predict the total pressure 
loss through the passage region. Since no external work is done on the flow, only 
viscous work, no streamline should leave the duct with a total pressure larger than 
the maximum inlet total pressure. Thus a total pressure loss coefficient referenced 
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to the maximum inlet total pressure must always be greater than zero. This 
calculation for the exit flow in the present case is shown in Fig. 7. Significant 
regions of total pressure greater than the maximum inlet total pressure are shown as 
the crosshatched region with the worse case occurring near the hub. The large 
increase in total pressure is apparently due to the pressure decomposition and 
approximation procedures used in the present approach. Since the dynamic head 
was very small for the current test case (and thus the dynamic pressure q was quite 
small) it is also clear that the numerical truncation errors of the current approach 
could be contributing to this anomalous behavior. A subsequent calculation of 
higher Mach number (and thus higher q r ) will be discussed later to give more 
insight into this problem. 

A final consistency argument can be made with regard to the transverse 
pressure gradients obtained from the current approximate analysis. Since the two 
transverse momentum equations are not used in the analysis; they may be used as 
"a postori" tests of the accuracy of the solution. On the endwall, the transverse 
momentum equation reduces to 


1 _ 3P 0 y as 

h 2 dy 2 h 3 ay 3 (4.1) 

\ 

thus relating transverse pressure gradient to vorticity gradient. Eq. (4.1) has 
been integrated along the endwall using the value of £ from the solution of 
equation (3.3a). The resulting pressure distribution is compared_with that 
obtained from the present approximate model, Eq. (10.2.16), in Figure 8 . Clearly 
the transverse pressure field obtained from Eq. (10.2 . 16) is not in balance with 
the transverse momentum equation and additional work may be required on this 
approximate model of the flow field. 

Large Inlet Dynamic Pressure 

The baseline test case study for laminar flow through the turbine duct 
presented above was performed for an inlet Mach number M = .1. Therefore the 

dynamic pressure at the inlet was less than one percent of the static pressure so 

that small errors in calculating static pressure could lead directly to significant 
errors in the static pressure and total pressure coefficient. Therefore in order 
to assess the level of accuracy being achieved in predicting pressure levels, the 
baseline test case was rerun with a higher inlet Mach number (M = 0.25) to 

produce an inlet dynamic pressure 7 times that of the base case (Case 2 Table I) . 

The calculated results for total pressure loss distribution are shown on Fig. 9. 
Again regions of significant gain are encountered with only marginal improvement 
over the lower Mach number results. It is apparent that the total pressure loss 
through the duct is not yet adequately represented by this approximate flow field 
model. 

A final note should be made concerning this case. Since the vorticity trans- 
port equation (3.3) is incompressible and the complete energy equation is not used 
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(Eq. 3.12), the calculation is more nearly an incompressible flow calculation 
although density and temperature variations occur. Thus when the inlet Mach 
number was raised to .25 and subsequently to .30, the calculation failed half 
way through the duct at approximately the minimum area ratio. The indicated Mach 
numbers at this last station were M = .87 for the Mr = .25 case and M = 1.4 for the 
M r = .30 case. These observations indicate that calculations at higher Mach 
numbers where compressibility is significant may not yet be reliable. 
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4.2 Circular Arc Duct 


The case of flow in a circular arc duct represents an important problem 
because for some configurations comparisons can be made with experimental data 
published by Mori (Ref. 14). The general geometry of this case is shown in Fig. 

10 along with the definition of the principal parameters effecting the solution. 

Of these parameters, the most important is Dean’s number which represents the 
ratio of the viscous force to transverse static pressure force. The flow para- 
meters for these cases, designated 3 and 4 in Table I, represent the flow in a 
square duct with different step size so that the flow turns through 90° and 220°, 
rspectively. Of these parameters, it should be noted that Re<j is based on 
hydraulic dyameter and R e on gap width. 

The solution for the developing profile of the streamwise velocity at the 
plane of symmetry for case 3 is shown on Fig. 11. For this case, as in previous 
cases, the integration used a 20 x 20 mesh in the transverse plane and 24 stream- 
wise stations for a streamwise steps size A0 = 3.68 . The growth of the wall 
boundary layer is clearly evident and the free stream shows a pronounced velocity 
gradient with a peak on the suction side of the duct. This velocity gradient is 
much larger than would be obtained by simple radial equilibrium of two dimensional 
flow, thus indicating substantial secondary flow effects. 

A second case was run for the circular arc duct (case 4) which was the same 

as case 3, but with a larger step size A0 = 11.05 . For this case the flow made 

a turn of 220 . The streamwise velocity distribution at the plane of symmetry at 
220 is shown at Fig . 12 where it is compared with the experimental data of Mori 
(Ref . 14) . This comparison is quite good in the mid-stream region and indicates 
that the essential secondary flow phenomena is well represented by the vorticity 
transport equation. Considering the assumptions made for the vorticity transport 
equation (Eq. 3.3), it may be concluded that a large region of the secondary flow 
in the midspan region is essentially inviscid and can be represented fairly well 

by Eq. 3.7. Near the wall, however, insufficient data exists to verify the 

approximate model used in these calculations. In particular, as shown on Fig. 

12, a noticable difference exists between the calculated velocity and measured 
velocity near the pressure side of the duct. This difference may be due to the 
inviscid boundary conditions on the vorticity transport equation. 
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5.0 CONCLUSIONS 


Based on the results presented in this study it can be concluded that : 

1. - The explicit use of the vorticity transport equation to describe 
secondary flows and its interaction with the primary streamwise flow leads to 
the appearance of a clearly defined passage vortex. 

2. Approximations used in the analysis, including pressure splitting and 
inviscid like boundary conditions on the vorticity transport equation may effect 
quantitative predictions . 

3. Insufficient experimental data exists to verify the approximate flow 
model in regions of large secondary flow velocities . 

In particular, the accuracy of the flow conditions at and near the wall 
would be expected to be sensitive to the approximations made in constructing 
and solving the local flow model. 
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7.0 LIST OF SYMBOLS 


C 

Cp, 


Cu 


D 

E 

G 

^ 1*^2 9 ^3 
K 


M 

P 

R 


*S 

Re 

S 

T 

UjVjW 

X,Y,Z 


Speed of Sound 
Specific Heats 
Pressure Coefficient 
Hydraulic Dynameter 
Internal Energy CpT 
Gap 

Metni Scale Coefficients 
Dean’s Number 
Reference Length 
Mach Number 
Static Pressure 
Radius of Curvature 
Gas Constant 
Reynolds Number 
Span 

Temperature 
Velocity Components 
Coordinates 


Y 

K 


h 

P 


♦a 

'i's 


Ratio Specific Heats 
Streamwise Vorticity 
Principle Normal Vorticity 
Density 

Secondary Velocity Potential 
Secondary Stream Function 


Subscripts 


r Reference Condition 

I Inviscid Flow 

S Secondary Flow 

V Viscous Flow 
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8.0 FIGURES AND TABLES 
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TABLE I TEST CASE CONDITIONS 



Case No. = 

1 

2 

3 

4 

Reynolds Number R e 

= 

1000 

1000 

940 

940 

Mach Number M 

= 

.1 

.25 

.1 

.1 

Temperature °R T r 

= 

416. 

416. 

416. 

416 . 

Pressure atm P r 

= 

.25 

.25 

.25 

.25 

Density slug/ft^ p r 

= 

7.66-04 

7.66-04 

7.66-04 

7.66-04 

Viscosity ft^/sec V r 

= 

.0165 

.0143 

.0165 

.0165 

Velocity ft/sec U r 

= 

100. 

250. 

100. 

100. 

Prandtl No. P^ 

= 

1.0 

1.0 

1.0" 

1.0 

Span/Gap S/G 

= 

3.0 

3.0 

1.0 

1.0 

R/D 

= 



14.0 

14.0 

re d ... 

= 

- - 

- - 

9401 ‘ 

940. 

Dean's Number K 

= 



251. 

251. 

ARC Length (deg) 0 

= 



90. 

220. 
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FIG. 3 


STREAMWISE VORTICITY DISTRIBUTION AT (TURBINE) DUCT EXIT 

PLANE OF SYMETRY 


2 Z/S 

4- A 0 014 
6- 0 0 032 
8- 0 0.063 
10-0 0.113 
12 - □ 0 194 
14-0 0.318 
16-/3 0 497 
18 - V 0.731 



ORIGINAL PAGE It 
OP POOR QUALITY 



0 0.2 0 4 0 6 0 8 1.0 


GAPWISE DISTANCE - Y/G 

78 - 04 - 187-4 


8-5 





SPANWISE DISTANCE Z/(S/2) 




ORIGINAL PAGE IS 
OP POOR QUALITY 


FIG. 5 


SPANWISE VELOCITY DISTRIBUTION FROM PRESSURE 
TO SUCTION SIDE AT (TURBINE) DUCT EXIT 




GAPWISE DISTANCE Y/G 
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FIG. 6 
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LOSS DISTRIBUTION AT (TURBINE) DUCT EXIT 
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DEVELOPING PROFILE IN CIRCULAR ARC DUCT 
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FIG. 12 
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FOR CIRCULAR ARC DUCT 
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9.0 PROGRAM DESCRIPTION - 


The Turbine End-Wall Computer Program (TEWCP) consists of two separate computer 
codes: 1) coordinate calculation (GDC00R) , and 2) three-dimensional viscous calculation 
(PEPSIP) . The GDC00R code calculates an orthogonal coordinate system based on 
conformal mapping (Ref. 6). In addition, a separate code (GCC00R) is provided 
to generate coordinates for circular arc ducts. The PEPSIP code is described in 
this report and the coordinate calculation in Ref. 6. 


9 . 1 Run S tream 

The PEPSIP computer code has no data card input. Each case must be treated 
separately and a new program constructed for each case. Most input parameters 
are set in subroutine INPUTS. Therefore the runstream for UNIVAC computers, requires 
compilation and construction of a new absolute code. Program GDCC^0R and GCC00R 
have data card input. The PEPSIP code also uses input data generated by GCCCfl^R 
describing the coordinates and stored on temporary or permanent files. Also 
depending on the PEPSIP options, there is a start/restart data file and an output 
plot data file. Therefore the runstream depends on the selected options. 

Three runstreams are described below. Run 1 shows the normal runstream where 
the general coordinates and viscous calculation are completed in the same run. 

Run 2 shows the runstream for circular arc ducts. Run 3 shows the runstream in 
which the coordinates have been previously calculated and stored on a permanent 
file. 
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Run 1 General Coordinates 


@ASG, A QUAL*PL0TD. , F111800 
OUSE 8 , QUAL *PL0TD . 

0ASG,T 9,0/ 50000/TRK/50000 
@ASG,T 10. 

@ASG,A QUAL*RESTART. ,D/100000/TRK/150000 
OUSE 11 , QUAL*RESTART . 

@ASG,T 12 ,D/5000/TRK/5000 
@ASG,T 13 ,D/ 5000/TRK/5000 
@ASG,T 14 ,D/ 5000/TRK/5000 

@ASG,A QUAL*ADD. 

0 COPY, A QUAL *ADD . GDC 0R , TP F $ . 

OXQT GDC0R 


\ input data ( 

Jnote: NDRUM = 10 | 

@ASG,A QUAL*PEFSIP. 

@COPY QUAL*PEPSIP,TPF 
@F0R,S INPUTS , INPUTS , INPUTS 

C make temporary changes to 
\ subroutine INPUTS 

@MAP , SEN MAPB,ABSA 
@XQT ABSA 

OFIN 


Plot data file 


Restart file 


Setup Coordinate 
Calculation 


Setup Viscous 
Flow. Calculation - 
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Run 2 Circular Arc Duct 


@ASG, A 
@USE 

QUAL*PL0TD. , F///800 
8 , QUAL *PL0TD . 

@AS G , T 
@AS G , T 

9, D/50000/TRK/50000 

10. 

@ASG,A 

@USE 

QUAL*RESTART. ,D/100000/TRK/150000 
11 , QUAL* RE START . 

@ASG, T 
@ASG,T 
@ASG,T 

12 s D/5000/TRK/5000 
13 ,D/ 5000/TRK/5000 
14 ,D/5000/TRK/5000 

@ASG, A 
0 COPY, A 
@XQT 

QUAL* ADD. 

QUAL* ADD. GDC0R,TPF$ . 
GDC0R 

\ input data / 

Jnote: NDRUM = 10 | < 

{?ASG,A 

@COPY 

@F0R,S 

QUAL*PEPSIP. 
QUAL*PEPSIP ,TPFS . 
INPUTS , INPUTS , INPUTS 


l make temporary changes to 
( subroutine INPUTS 

@ MAP, SEN MAPB,ABSA 
@XQT ABSA 

(§ FIN 


Plot data file 


Restart file 


Setup Coordinate 
Calculation 


Setup Viscous 
Flow Calculation 
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Run 3 Stored Coordinates 


@asg,a 

@USE 

QUAL*PL0TD.,F///8OO 
8 , QUAL*PL0TD . 


Plot data file 

@ASG,T 

9 jD/50000/TRK/ 50000 



@ASG,A 

@USE 

QUAL*C00RE. 

10 ,QUAL*C00RF 

1 

Coordinate data file 

@ASG, A 
@USE 

QUAL*RESTART . ,D/100000/TRK/150000 
11 , QUAL*RES TART 



@ASG,T 

@ASG,T 

@ASG,T 

12 , D/5000/TRK/5000 

13 , D/ 5000/TRK/ 5000 

14 , D/ 5000/TRK/5000 



@ASG, A 
@C0PY 
@F0R, S 

QUAL*PEPSIP. 
QUAL*PEPSIP. , TPF$ 
INPUTS , INPUTS , INPUTS 


Setup Viscous 
Flow Calculation 


make temporary changes to 
subroutine INPUTS 


@MAP,SEN MAPB.ABSA 
@XQT ABSA 

@FIN 
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9.2 Input for General Coordinates (GDC00R) 


Card No. 1 Title Card 


Name 

Col. 

Format 

Comments : 

Title 

1-72 

12A6 

Any Alpha Numeric Characters 

Card No. 2 

Option Card 



Name 

Col. 

Format 

Comments : 

JL(NS) 

1-2 

12 

No. Streamwise Stations 

KL 

3-4 

12 

No. Streamlines 

KDS 

5-6 

12 

= 1 

JLPTS 

7-8 

12 

No. Input Data Points 

IDBG6 

9-10 

12 

= 0 No Print, = 1 print coordinates 

NDRUM 

11-12 

12 

Disc Storage Unit No. =10 

LFILE 

13-14 

12 

= 0 Tape Storage 

DDS 

15-24 

F10.0 

Mesh Distortion Parameter 

DUCTI (1) 

25-34 

F10.0 

Duct Length (ft) 

DUCTI (2) 

35-44 

F10.0 

No. Knots 

Card No. 3 

through 2 JLPTS /8 

+ 5 


Name 

Col. 

Format 

Comments : 

RD1I (J) 

1-80 

8F10.0 

0D Duct Radius (ft) (J=l, JLPTS) 

RD2I (J) 

1-80 

8F10.0 

ID Duct Radius (f t) (J=l, JLPTS) 


9.3 

Output for General Coordinates 

USCALE 

u s 

Metric Normalizing Factor 

EPS (1) 

e l 

Roberts Transformation 

Parameter 
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9.4 Input for Circular Arc Duct 


Card No. 1 Title Card 


Name 

Col. 

Format 

Comments : 

Title 

1-72 

12A6 

Any Alpha Numeric Character 

Card No. 2 

Option Card 



Name 

Col. 

Format 

Comments : 

JL 

1-10 

110 

No. Streamwise Stations 

KL 

11-20 

110 

No. Streamlines 

KDS 

21-30 

110 

= 1 

JDBG6 

31-40 

110 

= 0 No Print, = 1 Print Coordinate 

NDRUM 

41-50 

110 

= 10 

DDS 

51-60 

F10.0 

Mesh Distortion Parameters 

Card No . 3 

Parameter Card 



Name 

Col. 

Format 

Comments : 

DUCTI (1) 

1-10 

F10.0 

d l 

DUCT1 (2) 

11-20 

F10.0 

d 2 

DUCTI (3) 

21-30 

F10.-0-- 

d 3 

DUCTI (4) 

31-40 

F10.0 

d 4 

DUCTI (5) 

41-50 

F10.0 

d 5 

DUCTI (6) 

51-60 

F10.0 

d 6 

DUCTI (7) 

61-70 

F10.0 

d 7 

DUCTI (8) 

71-80 

F10.0 

d 8 


9.5 

Output for Circular Arc 

Duct 

USCALE 

Us 

Metric Normalizing Factor 

EPS (1) 

£ 1 

Roberts Transformation 

Parameter 
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9.5 Input 

for PEPSIP 

Name 

Symbol 

Comment 

CP 

C p 

Specific Heat (f tlb/slug/°R) 

EPS (1) 

e l 

Roberts Parameter Gapwise 

EPS (2) 

£ 2 

Roberts Parameter Spanwise 

EZERO 

Eo 

Reference Energy lAIr-l)!^ ) + 1/2 

GAMMA 

Y 

Ratio Specific Heats 

IRS TIN 

= 0 

No Restart 


> 0 

Starting Station 

NFILE 

= 1 

No Restart 


> 0 

Sequence Number from Last Printout 

NS 


Last Station 

NSAVE 

- 0 

No Restart 


> 0 

Sequence Number from Last Printout 

PRL 

Erl 


PRT 

P RT 

jNote Set Equal to One in Code 
Gas Constant (f tlb/slug/°R) 

RG 

R 

RZERO 

Pr 

Reference Density 

S0UND 

C 

Reference Speed of Sound (ft/sec) 

TSTAG 

T s 

Reference Stagnation Temperature (° 

WALL 

T 

x w 

Wall Temperature (°R) 

TZERO 

T r 

Reference Temperature (°R) 

USCALE 

u s 

Metric Normalizing Factor 

UZERO 

u r 

Reference Velocity 

VISC0S 


Reference Kinematic Viscosity 

YS (2,1) -YS (1,1) 

g 

Inlet Gap (ft) 

YS (2,2) -YS (1,2) 

S/2 

Inlet Span (ft) 

X(2)-X (1) 

AX 

Initial Step Size 

YZERO 

Y r 

Reference Length (Set = 1.0) 


9.6 Output For PEPSIP 

Name 

Symbol 

Comment 

X 

X 

Streamwise Station 

Y 

Y 

Gapwise Station 

Z 

2 

Spanwise Station 

U-VEL 

u/u r 

Streamwise Velocity 

VC0MP 

v/u r 

Gapwise Velocity 

WC0MP 

w/u r 

Spanwise Velocity 

RH0 

p/p r 

Density 

V0R 

Streamwise Vorticity 

STREAM 

Y/(V r U r ) 

Secondary Stream Function 

CP/ 2 

(P-P r /(l/2P r U r 

^Pressure Coefficient 

TSTAT 

T/T r 

Static Temperature 
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9.7 PEPSIP/GDC00R Interface 


Let (X,Y,Z) be the generalized orthogonal coordinates used in the PEPSIP 
code and (s,n) be the two dimensional generalized orthogonal coordinates computed 
by GDC00R as described in Ref. 6. The coordinates are matched as follows: 


Y r h x d x = r r dS/V (9.7.1) 

Yj-hjjdy = r r dn/V (9.7.2) 


where r^ is the normalizing length scale used on GOC00R and V is the inverse of 
the metric scale coefficient, Y r is the length scale in PEPSIP. Then (x;s) is the 
streamwise coordinate matching pair and (y;n) is the gapwise coordinate matching 
pair. The spanwise coordinate is generated by extension 

Y r h z d z = dt/V(s,n) (9.7.3) 

The matching of length scales and metric coefficients is done on the coordinates 
defined at the inlet. Thus with V a metric normalizing factor, duct arc lengths 
at the inlet are given by, 

(9.7.4) 

r r 7-V r - - - ''(9.7.5)' 


(9.7.6) 


X(2)-X(l) = Y^AX = rj.AS/Vj 


n-p 


YS(2,1)-YS(1,1) = Y r J h y d y = r r 

y± 


/ 

n l 


dn 


YS(2,2)-YS(1,2) = Y r j h ^ f- 


/ 


Then we have using (9.7.2) and (9.7.5) with V r constant, 

- YS(2,1)-YS(1,1) . ( V r ) 

K~ h z ~ V Y~ ~ ~ y \ T~ I 


(9.7.7) 


Finally since V is the potential flow velocity, the incompressible pressure 
coefficient is 

/ \2 
V 


C = 1-1- 


(9.7.8) 
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Then with the PEPStP length scale Y r = 1, the coordinates are uniquely 
defined by 

X(2)-X(l) = first step size in marching direction (ft) 

YS(2,1)-YS(1,1) = inlet gap (ft) 

YS(2,2) - YS(1,2) = inlet height (ft) 

U = V = Metric normalizing factor. 

IT 
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10.0 APPENDIX 


An Approximate Analysis for Three-Dimensional Viscous 
Subsonic Flows with Large Secondary Velocities 

by 

W. R. Briley and H. McDonald 


The following section contains results from an analytical investigation of 
turbine endwall flow and heat transfer, performed by Scientific Research 
Associates, Inc., under subcontract to United Technologies Research Center. 
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FOREWORD 


The following report contains results from an analytical investigation of 
turbine endwall flow and heat transfer, performed by Scientific Research Associates, 
Inc., under subcontract to United Technologies Research Center. The investigation 
was initiated while the authors were affiliated with United Technologies Research 
Center and completed under Subcontract. 
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An Approximate Analysis for Three-Dimensional Viscous 
Subsonic Flows with Large Secondary Velocities 

INTRODUCTION 

An approximate analysis is developed for a class of three-dimensional subsonic 
viscous flows at high Reynolds number, which have a predominant primary flow 
direction with secondary flow in transverse planes. A key restriction of the 
analysis is that once a primary flow direction is established, reversed flow in the 
primary flow direction (flow separation) does not occur. As a consequence, the 
analysis is not intended for application in geometries having sharp curvature in the 
primary flow direction. However, the analysis is specifically designed for flow 
configurations having large turning as a result of strong geometric curvature, 

R = 0(1) but not sharp R = 0(e) curvature in the primary flow direction, R being 
the geometric radius of curvature. Such flows are common in turbomachinery applica- 
tions, where internal flow passages are usually aerodynamically smooth. The 
problem of primary interest here is the flow in turbine blade passages, and 
particularly the strong secondary flows and corner vortices associated with the 
endwall regions. 

The present analysis synthesizes concepts from potential flow theory, secondary 
flow theory, and from an extension of three-dimensional boundary layer theory in a 
manner which allows efficient solution by forward-marching numerical integration 
techniques, given only knowledge of the potential flow. Repeated and costly 
(iterative) sweeping of the entire three-dimensional flow field is not necessary 
to obtain a valid approximation to the flow. As a consequence, numerical solution 
of the derived governing equations will generally require considerably less comnuter 
time than would a solution of the full three-dimensional Navier-Stokes equations. 

An £_ priori potential flow solution nevertheless accounts at least anproximately 
for the important elliptic influence of downstream boundary conditions and 
determines the primary flow direction. One might consider a possible variation 
of the analysis which would employ an ja priori invlscid flow not from a potential 
flow solution but from a solution of the rotational inviscid flow equation. This 
would result in a considerable increase in computational effort since for the same 
number of mesh points solving the rotational inviscid flow equations takes only 
slightly less effort than solving the full Navier-Stokes equations. On the other 
hand, much less effort is required to solve the potential equation. Finally, the 


10-3 



analysis recognizes that a pressure field determined solely from potential flow 
theory is probably an inadequate representation of the pressure field when large 
confined secondary flows are present. 
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BACKGROUND 


Because of their complexity, and particularly the interaction which occurs 
between primary and secondary flows and between viscous and inviscid regions, three- 
dimensional viscous flows, particularly confined flows m highly curved ducts, have 
been extremely difficult to analyze, short of numerical solution of the comnlete 
Navier-Stokes equations. This latter alternative is of course definitive but, in 
three dimensions, it requires a great deal of computational effort. The present 
analysis, although approximate, is much more economical in terms of computer time. 
Approximate methods, wrhich treat the flow partly by forward marching and partly as 
fully elliptic, have been proposed recently by Dodge [1] and by Pratap & Spalding [2]. 
These methods require iterative sweeping of the entire three-dimensional flow 
field, however, and consequently the computer time required tends to be comparable 
with that needed for solution of the complete Navier-Stokes equations, without 
approximation. The present analysis is an outgrowth of and to some extent a 
synthesis of previous work involving widely differing approaches. Highlights from 
this previous work are briefly reviewed here. 

Rotational Inviscid Flow Theory 

Rotational inviscid flow” theory has provided considerable insight into the 
behavior of some secondary flows, and in many cases can provide quantitative predic- 
tions. For flows w 7 ith large turning. Squire & Winter [3] first proposed a wridely 
used inviscid "secondary flow approximation" which governs the generation of streamwise 
vorticity due to turning. Secondary flow theory and its applications have been 
reviewed by Hawthorne [4], Horlock & Lakshminarayana f5] and Lakshminaravana & 

Horlock [6]. Stuart & Hetherington [7] have computed solutions to the full incompres- 
sible, rotational inviscid equations of motion for flow through a curved duct. 

Treatment of viscous effects, however, has been difficult, although the relevant 
governing equations are known. Hawthorne [4] points out the necessity of defining 
a "cut-off velocity" near walls, to exclude the viscous region, and that this can 
have a major effect on the predictions of secondary flow theory. Techmaues for 
computing three-dimensional boundary layers are available (cf, Nash & Patel [8]); 
however, three-dimensional boundary layer theory is not valid in viscous corner 
regions, and the patching or interfacing of boundary layer and rotational inviscid 
regions is not at all straightforward. 
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Extended Boundary Layer Analysis 

Patankar & Spalding [9], Caretto, Curr & Spalding [10], and Briley [11] 
have suggested numerical methods for solving approximate governing equations which 
represent an extension of concepts from three-dimensional boundary layer theory. 

In these studies, solutions were computed for laminar incompressible flow in straight 
ducts with rectangular cross sections. The governing equations were solved by 
forward-marching integration in a primary flow coordinate direction while retaining 
viscous stresses in both transverse coordinate directions as opposed to only one 
direction for three-dimensional boundary layer theory. Viscous diffusion in the 
primary or streamwise direction is neglected as in conventional three-dimensional 
boundary layer theory. Unlike conventional boundary layer theory, however, the 
approximate governing equations are solved in the inviscid region as well as the 
viscous region. In addition to neglecting streamwise diffusion, the pressure is 
divided into an inviscid pressure "imposed", for example, by an irrotational flow, 
and a pressure correction due to viscous effects. The "viscous" pressure is treated 
differently in the streamwise and transverse momentum equations (an inconsistency) 
to permit solution by forward marching integration. Despite the inconsistent 
treatment of pressure gradients, these so-called "parabolic flow" methods gave 
reasonable predictions for flow in straight rectangular ducts. Ghia, Ghia & 

Studerus [12'] also obtained reasonable predictions for flow in a "polar" duct whose 
cross sections are annular sectors. Recently, as part of a comprehensive studv of 
entry flow m straight rectangular ducts, Rubin & Khosla [13] presented their own 
analysis and numerical method for the "fully viscous region". Their method is 
an extended boundary layer method in the present terminology. 

Although the general approach of [9-13] produces reasonable results for flow 
in straight ducts, only very small secondary flows are present in such cases. 
Nevertheless, this approach has been used to compute moderate secondary flows 
generated by a moving wall [9], by a transverse buoyancy force [11], and in mildly 
curved pipes [14] and curved rectangular ducts [15], In the fully-developed 
flow region, these results generally agree with other fully-developed solutions 
computed without approximation. Very little is known about the accuracy of these 
solutions in the developing flow region. It is not expected that these analyses 
as they presently exist would produce reasonable results for developing flows in 
ducts with la'rge turning .producing strong secondary flows, the problem of present 
concern. As a consequence, the present investigation was undertaken and had as its 
main goal the synthesis of results from "extended" boundary layer theory and classical 
secondary flow theory. 
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ANALYSIS 


Outline 

The present analysis is based on the decomposition of the overall velocity 
vector field U into "inviscid" components U^. and U g which satisfy slip conditions, 
and a "viscous" component which corrects for no-slip conditions at, wall boundaries. 
The overall velocity is determined from the superposition 

U=D 1 +U s +U v (10.2.1) 


Elements of the velocity decomposition are briefly described below: 

(1) Uj is the irrotational velocity field determined from a potential flow solution 
in the geometry under consideration. Depending on the particular application, 
approximate downstream conditions such as the commonly used "Kutta condition" 

may be required. The direction of U is then taken as the "primary flow direction". 
Curvilinear' surfaces of constant potential may be taken as the transverse 
"secondary flow" planes, although this is not essential. 

(2) U is a secondary flow velocity field which is further decomposed into a 
solenoidal component and an irrotational component U^. by definition has 
zero divergence and is determined from an approximate application of secondary flow 
theory, under assumptions less restrictive than those of the now classical 
"secondary flow approximation" of Squire & Winter [3]. Briefly, is obtained 

from the solution of an approximate equation governing the growth and development 
of streamwise vorticity, followed by approximate vector potential calculations in 
curvilinear transverse secondary flow planes. The troublesome "velocity cut-off 


problem" of applications of inviscid secondary flow theory, mentioned by Hawthorne 
[4], is avoided by the retention of viscous terms. The irrotational component 

i 

ensures that the overall velocity field satisfies continuity and is obtained 
from approximate scalar potential calculations in curvilinear transverse secondary 
flow planes. 

(3) U is a viscous velocity correction for no-slip conditions, computed under the 

influence of "imposed" pressure gradients determined by the superposition of 

and U . 
s 

Computationally, U^. is determined by an a_ priori elliptic potential flow 
calculation, whereas U g and are computed simultaneously by a coupled forward- 
marching numerical integration procedure. It is unnecessary to sweep the three- 
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dimensional flow field iteratively to obtain a valid approximate solution, although 
such a procedure is possible and may possibly improve the flow predictions. A key 
feature of the present analysis is that the pressure gradients associated with 
large secondary flows are determined a posteriori using approximations to these 
secondary velocities. This is necessary since although velocity is sufficient 
to determine pressure gradients from the momentum equations, a knowledge of 
pressure is insufficient to determine velocity, since velocity must satisfy 
continuity in addition to the momentum equations. Finally, although the present 
analysis can be easily extended to give a more comprehensive treatment of compres- 
sibility effects, at present compressibility is accounted for only through the 
use of the perfect gas equation of state in the imposition of the streamwise 
pressure gradients. 

Primary Flow and Transverse Secondary Flow Planes 

The primary flow is taken as the potential flow past the geometry m question. 

Thus, 


°i ‘ 

and from continuity, 

V-0 X = V 2 <£ r = 0 


( 10 . 2 . 2 ) 


( 10 . 2 . 3 ) 


The potential flow U^. can be determined either analytically or by numerical 
computation in some convenient coordinate system. The direction of is taken as 
the primary flow direction. The unit vector i in the primary flow direction is 
given 'by 


'p- U I /U ! 


( 10 . 2 . 4 ) 


where velocity magnitude is denoted U and defined by 


U = 



(D-D ) l/2 


( 10 . 2 . 5 ) 


It is also necessary to specify transverse secondary flow planes, which may be 
curvilinear. The unit vector normal to the secondary flow planes is denoted i^. 
Curvilinear planes of constant potential <f> x are an obvious but not essential 


10-8 



ORIGINAL PAGE IS 
OF POOR QUALITY 


choice for the secondary flow planes, and in this case, i 



Secondary Flow Velocity Components 


The solenoidal component of the secondary velocity U g is determined from 
an approximate application of secondary flow theory. From Horlock & Lakshminarayana 
[5] the equation governing the growth of streamwise vorticity for incompressible 
flow is 


ds \ v / vR v 2 l J 


( 10 . 2 . 6 ) 


where s is distance along a streamline, s is the unit vector tangent to a 
streamline, £2 is the vorticity vector, £ is vorticity in the direction of s, 

V is velocity magnitude, R is the local radius of curvature of a streamline, 
n is transverse vorticity in the direction of R, and v is kinematic viscosity. 
Equation (10.2.6) is exact and involves no approximation. 

To obtain E, and subsequently U^, (10.2.6) is solved subject to the following 
approximations : 

(1) C is taken to be vorticity in the direction of the primary potential flow 

streamlines, i.e., £ = £(i • s) 4 £ 

^ P 

(2) s is taken to be distance in the direction of + U s , i.e., 

a . (Qi+o s ) 7 (10.2.7) 

as = IDj + DJ 


Note that the inclusion of U in the secondary flow convective operator is a key 

s 

feature of the present analysis, since this is the mechanism whereby streamwise 
vorticity in a curved blade passage is transported into the suction surface - 
endwall corner where it becomes the well-known corner vortex. 

(3) V is taken to be the computed streamwise velocity magnitude, i.e., 

V = i • (U, + U + U ) . 
p i s v 

(4) R is taken to be the local principal radius of curvature of a primary flow 
streamline. 

(5) p is taken to be the transverse vorticity in the direction of R as computed 
from V. 

(6) The bracketed viscous diffusion term in (10.2.6) is replaced by an expression 
representing diffusion only in the transverse secondarv flow planes: diffusion 

in the streamwise direction normal to these planes is neglected. This is 
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accomplished by replacing the bracketed term in (10.2.6) by 


-T n -(VxVxf n i n ) = s * V 2 £2 


(10.2.8) S 


where C n is the projection of £ g ip onto i^. 


i.e. 


» ^s (l n 


V* 


The solenoidal velocity is determined from an approximate vector 
potential calculation based upon the classical definition of vector potential 
(e.g., Newell [16]): 


v = vx¥ , a- vxv = -v 2 ^ 


(10.2.9) ] 


Here, it is assumed that ¥ is perpendicular to the transverse secondary flow 
planes and thus has the same direction as i^. Thus, 

U^ £ Vx\// s T n (10.2.10)' 


Consequently, 


U, 


lies in the transverse secondary flow planes and 




= i n 1 t V X V X 


s'n^ 


( 10 . 2 . 11 )' 


The function is analogous to a stream function defined in the transverse planes 

at each streamwise or axial location. Equation (10.2.11) is elliptic in the transvers 

planes and contains no streamwise derivatives, as a result of the assumption 

Y = ib i . Note that ib and thus U. can be determined at each transverse plane 
s n s ip 

from a knowledge of secondary vorticity £ n in that plane, and no knowledge of 
at planes upstream or downstream is required, although itself depends upon 
upstream information. 

The irrotational component U. of the secondary flow velocity U is determined 

Y S , 

from the requirement that the overall velocity U, including the viscous correction 
Up., satisfies local continuity. Actually, U^. and satisfy continuitv by 
construction, however as will be explained, U v is determined from momentum 
considerations together with an integral continuity restraint, and consequently 
U v does not necessarily satisfy local continuity. To ensure local continuity, U^ 
is determined such that 


V-(U^+U v ) =o 


( 10 . 2 . 12 ) 
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Furthermore,, it is assumed that U lies in the transverse secondary flow nlanes 


and is derivable from a surface potential (j> g such that 


u 4> s7< ^s- i n (T n- V< # , s > 


(10.2.13) 


Consequently, <t> s is governed by 


7 S '#’s- 7 ' T iA' 7 i # , s>-- 7 - u * 


(10.2.14) 


which, given a knowledge of U^, is elliptic in the secondary flow planes. 


Viscous Velocity Correction 

The primary and secondary flow velocity components U T and I) satisfy slip 

X s 

conditions at solid boundaries. U corrects for the no-slip conditions, based on 

v 

solution of an approximate momentum equation under the influence of "imposed" 
pressure gradients. The momentum equation is given by 

(U V)U + — Vp = F(U) (10.2.15) 

P 

where F is the force due to viscous stress. Substituting decomposition formulas 
for velocity and pressure 


u = 0 I + u s +u v , P - Pj + p 5 + P v 


(10.2.16) 


there results 

[(D I + O s + D v )-v](D I + 0 s + u v ) + ~V( Pi + p s +p v ) = f( 0 i + u 3 + D v ) ( 10 . 2 . 17 ) 

Borrowing a concept from boundary layer theory, it is assumed that the pressure 
gradients to be imposed are determined by + U , which satisfy slip contiions. 

Thus, 

— V( p i + p s ) = F(U X +U S ) - [{Uj +U s )-v](U I + U s ) (10.2.18) 

It should be noted that although U satisfies slip conditions, U and the resulting 
pressure gradients are not really "mviscid", since viscous terms are retained 
in the streamwise vorticity equation, and the irrotational contribution depends 
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upon the viscous correction U^. Combining (10.2.17) and (10.2.18), an equation governing 

the viscous correction U is obtained: 

v 

[(r r + U s ) v]u v + [u v -v](u i +u s +u v ) + *^-Vp v = F<U V ) ( 10 . 2 . 19 ) 

i 

Further assumptions are necessary to permit solution of (10.2.19) by forward 
marching integration for subsonic flows. Streamwise diffusion in the direction 
normal to the secondary flow planes is neglected. In addition, for the primary 
flow direction, p^ is taken to be a function only of distance s^ in the nrimary 
flow direction. For external flows, p^ is taken as zero in the absence of elliptic 
viscous-inviscid interaction. For internal flows, p^ represents a mean viscous 
pressure drop determined by the requirement that the integral mass flow be conserved. 

The momentum equation for the primary flow direction becomes 

Tp- {[(Uj+u s )-v] u v + [u v -v](Uj + 0 s + u v ) + —Vp v (s p )- f'(u v )} = 0 (10.2.20) 

■where F 1 does not include streamwise diffusion. The integral mass flux condition 
for internal flows with no mass addition is 



= constant 


( 10 . 2 . 21 ) 


Assumptions for p^ appropriate for the transverse components of the momentum 
equation are the topic of a continuing investigation. However, a very simple 
(although highly approximate) method for determining the transverse components of 
U v can be devised from arguments based upon Johnston's [17] empirical observations 
of three-dimensional turbulent boundary layer velocity profiles having triangular 
polar diagrams. The triangular polar plots display a structure such that in the 
outer (inviscid) region, the skewing of the velocity is consistent with secondary 
flow theory, and in particular is related to the turning of the free stream in 
accordance with the Squire & Winter [3] secondary flow approximation. In the 
inner (viscous) region, the velocity profile is found to be collateral, such 
that v/u = constant. This suggests the following formula for computing the 
transverse components of U^: 


u v = (u l + 0 s ) 


ip-tu^Us+Uy) 

V(0 T +0J 


( 10 . 2 . 22 ) 
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The formula (10.2.22) has the following properties: (1) in the outer portion of 
the viscous region, where ip . TJ^ 0, the composite velocity U reduces to TT^ + U g , 
(2) the composite velocity U satisfies the no-slip condition and is collateral 
near the wall. Improved approximations for the near-wall region are possible 
but are not considered here. 


Energy Equation 

Application of the present theory to nonadiabatic flows requires solution of 
the energy equation, which can be written 

V-/3UE = V-crVT + 0-F + $> (10.2.23) 

Solution of (10.2.23) by forward marching integration requires only that F be 
replaced by F' and that terms representing the streamwise conduction of heat be 
neglected. 


Compressibility Effects 


The foregoing analysis has been formulated primarily for incompressible 
flows, for which the density p is constant. For relatively small Mach numbers, 
however, compressibility effects can be represented in an approximate sense by 
introducing the perfect gas equation of state p - pRT in the imposition of stream- 
wise pressure gradients. Eliminating temperature T using the temperature-enthalpy 
relation 


T = 



U-U \ 

2 j 


(10.2.24) 


the following relation for representing streamwise pressure gradients is obtained: 


i p .V( p i + p s +p v ) 



(10.2.25) 


which provides an auxiliary equation relating densitv, velocity and 

P'S 

A slight modification results if U • U is replaced by (i * U)^. A 
hensive treatment of compressibility effects is possible but is not 
here. 


total enthalpy, 
more compre- 
considered 
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APPLICATION TO A CURVED FLOW PASSAGE 


The general analysis of the previous sections is applied to the flow in a 
curved duct which is shaped like a turbine blade passage. A cross section of 
the duct geometry is shown in Fig. 1 and consists of two curved walls corresponding 
to the surfaces of turbine blades. The remaining boundaries of the flow passage 
are two parallel flat surfaces comprising the endwalls. In the particular 
case being considered, the geometry and inviscid flow are two-dimensional. 

However, the viscous flow is three-dimensional, particularly near the endwalls, 
where the transverse pressure gradients associated with the turning of the flow 
produce strong secondary flow toward the low pressure or suction surface side 
of the passage. A primary objective of this initial application is to demonstrate 
that the present analysis can predict endwall flows and the associated corner 
vortex region which are typical of that known to exist experimentally in 
geometries of this type. 

Coordinate System 

Curvilinear orthogonal coordinates x, y, z are constructed to fit the flow 

passage boundaries as shown in Fig. 1. Metric coefficients hi, h 2s h 2 are 

2 7 

defined such that incremental distance 6s is determined by (6s) = (l^dx)^ 

+ (h 2 <Sy-) 2 + (h^z)^. In planes parallel to the endwalls, orthogonal streamlines 

and velocity potential lines from a two-dimensional incompressible potential flow 

> 

analysis are utilized as the coordinate lines for constant v and x, respectivelv. 

The z direction can be regarded as Cartesian, and thus h^ = 1. The x coordinate 
is the primary flow or axial coordinate, and is associated with surfaces for 
which 4>j is constant. The y and z coordinates define transverse secondarv flow 
planes at any given x location. In this coordinate svstem, the potential flow 
direction ip and the normals i^ to the transverse planes coincide with i, the 
unit vector in the x direction. A two-dimensional incompressible potential 
flow analysis and computer program developed bv Anderson [18] was emnloved in 
the present investigation, without modification, to compute the necessary coordinate^ 
data and potential flow solution. 

Governing Equations 

Throughout the remaining discussion, all variables in the governing equations 

are nondimens ional , having been normalized by the following reference quantities: 

2 

distance, L r ; velocity, U r ; density, p r ; temperature, T r ; total enthalpy, U r ; 
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pressure, p^U^; viscosity, Here the subscript r denotes a reference quantity. 

This normalization leads to the following nondimensional parameters: Mach number, M; 

Reynolds number, Re; Prandtl number, Pr; and specific heat ratio, y. These 
parameters are defined by 


M = 


= U r /c , Re =p r U r L T /fj. r , Pr = c /i r /k , 7 = C / 


(10.3.1) 


where is the molecular viscosity, k is thermal conductivity, and c^ and c^ 

are the specific heats at constant pressure and volume. The reference speed of 

2 

sound, c, is defined by c = YR T , where R is the gas constant. 

g r’ g 

Since in the present application the coordinate lines coincide with the 
potential flow streamlines and potential surfaces, the following simplifications 
hold: Vj = Wj = u g = 0. The velocity decomposition (10.2.1) can thus be expressed 

as 


U = lu+jv + kw =i(u I + u v )+j(v s + v v ) + k(w s +w v ) 


(10.3.2) 


where i, j, k are unit vectors in the x, y, z coordinate directions, respectively. 
Under the stated assumptions, the streamwise vorticity equation can be written as 


U I I du \ v s h, / d£ du \ w s / d( du \ 


2 3h, ^ fj. / d 1 


h 2 dy U dz + Re \ ay h,h 

2 a y 

2 2 

+ v + w , and the vorticity is related to 
s S 

, _ , . 1 ( dw s 

dh 2 v s \ 

c n ' c “ h 2 \ dy 

<?z j 

velocities are given by 


' ! 


V s ” h 2 dy h ( 

dz 

s _ _[ 

ah|^ s 

w s " dz hjh 2 

dy 


+ h 


a 2 e 

iaz 2 


) 


(10.3.3a) 


(10.3.3b) 


(10.3.4a) 


(10.3.4b) 
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(10.3.5) 


where the scalar and vector potential functions (j> g and <|i are governed by 


I 


h 2 ay h h e ay ^ + h, az 2 ' ^ 


and 


d h d<f> 


s _ . 


ay h g dy 


+ h.h. 


-( 


• 2 az 2 

The streamwise momentum equation is approximated by 


dh g u v ah^y ah,h g w 
ax + dy + az 




a , 

' a i 

ah 

a 

* 

ax (h *P u 1 + 

ay + T7 

ay 

h^u(v s + v v ) + h,h 2 — 

pU ( W £ + w ) 

■ 




S V 

h 


a(p x +p s ) 

2 "ax" 


dp v (x) 


2 dx 


I 


a p. au h,h 2 


Re dy h 2 ay 


Re az 


a. 


au 

az 


The energy equation is approximated by 


•^(h 2 puE) + 

a 

ay 

■[V (v s +v v 

,)E] 

a 

+ aF 

- [- h , h i 

h, a , 

r H- 

i 

aE , 1 / 

' 

H- \ 

d r 

Re ay ' 

L Pr 


ay 2h 2 \ 

H- ~ 

R- / 

ay L 

, h . h s 

a 


_£L + _!_/ 



a r , 

Re 

az 

l p r 

az + 2 \ 


Pr / 

az L u 


The imposed pressure gradients are determined from 

a i ah 


h a(p i +p s } _ 2\ 

2 ax " ~ ax h 2 u i 


i 


ay h ( ay 


a „ ah 2 

h ,Vs- h 1 h 2 -ar u i w s + ''|-5r 


original 
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(10.3.6) 


(10.3.7) 


(10.3.8) 


(10.3.9) 



and related to the dependent variables using the following form of the auxiliary 
gas law relation (10.2.25) 


a 

ax 


2-i 

r 



dtPj+Pg) dp v (x) 
ax + dx 


(10.3.10) 


The integral mass flux relation is 



h £i pudydz 


= constant 


(10.3.11) 


The collateral flow assumption (10.2.22) leads to the following equations for 
computing w : 


V 


v 


w 


V 



(10.3.12a) 

(10.3.12b) 


Numerical Method 

The governing equations are replaced by finite-difference approximations. 

Three-point central difference formulas are used for all transverse spatial 

derivatives. An analytical coordinate transformation devised by Roberts [19] is 

employed as a means of introducing a nonuniform grid in each transverse coordinate 

direction, to concentrate grid points in the wall boundary layer regions. Second- 

order accuracy for the transverse directions is rigorously maintained. Two-point 

backward difference approximations are generally used for streamwise derivatives, 

although this is not essential. The streamwise vorticity equation (10.3.3) is 

decoupled from other equations in the system and linearized with respect to £ bv 

an ad^ hoc process consisting of lagging quantities not yet available at the 

implicit level. The resulting implicit difference equation for £ is solved using 

a scalar ADI scheme based on the technique of Douglas & Gunn [20] for generating 

ADI schemes as perturbations of fundamental implicit schemes. Equations 

(10.3.5 - 10.3.6) for the scalar and vector potential functions (j>_ and ib 

*=> s 

are elliptic in the transverse planes and are solved, given values for the 
right-hand sides, using scalar iterative API. Specifically, the Douglas- - 
Gunn perturbation of the Crank-Nicolson scheme is used. The streamwise 
momentum (10.3.7), energy (10.3.8) and auxiliary gas law relation 
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(10.3.10) are solved as a coupled system using linearized block ADI. 

A detailed discussion of the block ADI scheme used has been given by McDonald & 
Briley [21] and Briley & McDonald [22], The general approach is to linearize 
the implicit equations by formal expansion about the solution at the most recent 
axial location. Terms in the difference equations are then grouped by coordinate 
direction, and the Douglas-Gunn [20] technique is used to generate an ADI scheme. 
The resulting difference equations can be written in block-tridiagonal or a closely 
related matrix form and solved efficiently by block elimination techniques. 

A summary of the overall algorithm used to advance the solution a single 

axial step follows. It is assumed that the solution is known at the n level 

n 
x 

( 1 ) 

( 2 ) 


n , . , . , n+1 

x and is desired at x 


The streamwise vorticity equation (10.3.3) is solved using scalar ADI to obtain 
r n+l 


The vector potential equation (10.3.5) is solved using scalar iterative ADI 
to obtain ij; n+ \ 


( 3 ) 


n 


( 4 ) 


Values for v and w are computed using equation (10.3.4). Values of <j> 

5 S 3 

obtained from the previous step are used. Values for v and w^ are 

computed from equations (10.3.12). The imposed pressure gradients 

3(p ] - + p g )/3x are computed from equation (10.3.9). 

A value for the mean viscous pressure drop dp^(x)/dx is assumed. Initially, 
the value from the previous step is used. 


(5) Equations (10.3.7), (10.3.8), & (10.3.10) are solved as a coupled system using 


block ADI to obtain values of u 


n+1 n+1 


and E 


,n+l 


In general, the 


( 6 ) 


integral mass flux relation (10.3.11) will not be satisfied. 

Return to step (4) and repeat this process ' iteratively using the standard 
secant method [23] to find the value of the mean viscous pressure drop which 
leads to u n+ ^ and p n+ ^ which satisfies the integral mass flux relation (10.3.11). 
The secant method was found in practice to converge to five figures on the 
third iteration. 

(7) Using values now available for 3u/3x, the scalar potential equation (10.3.6) 

is solved using scalar iterative ADI to obtain <}> n+ ^ for use in the next 

3 

axial step. 

To start this procedure, values for <j> n are needed. These are obtained by 
solving steps (4-7) above, under the assumption that the transverse velocity 
components are zero. 
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Computed Solution 

A sample laminar flow calculation was performed, to evaluate the overall 
method in its present state of development and thus explore the potential 
of the present method for making detailed predictions of the flow and heat 
transfer in turbine blade passages. The solution was computed for the geometry 
of Fig. 1, and the Reynolds number based on the duct width at the inlet was 1000. 
The ratio of the distance between endwalls to the distance between blade surfaces 
at the passage entrance was 3.08. Blade boundary layers at the upstream starting 
plane were assumed to have a thickness of 0.05 times the distance between blade 
surfaces. The endwall boundary layer was taken to be twice this thickness. 
Adiabatic walls were assumed, and the following flow parameters were used: 

Mach number M = 0.1, Prandtl number Pr = 1.0, specific heat ratio y = 1.4. 

A 20 x 20 grid was used for the cross section, and 24 axial steps x?ere computed. 

The ability of the present analysis to predict the formation and development 
of a corner vortex is of primary interest here. The existence of a corner vortex 
for this type of flow problem has been established for some time (e.g., Herzig, 
Hansen & Costello [24]). The flow structure consists of a strong secondary flow 
near the endwall and toward the low pressure or suction surface, caused by 
transverse pressure gradients associated with the turning of the primary flow. 

Upon reaching the suction surface, the secondary flow proceeds up the suction 
surface a short distance and then back into the primary flow, forming a stream- 
wise vortex m the corner formed by the intersection of the suction surface and 
the endwall. Recently, more detailed experimental measurements and observations 
have been reported by Langston, Nice & Hooper [25] and Sjolander [26]. The 
authors are unaware of any previous numerical calculation w T hich correctly predicts 
the formation and subsequent development of the corner vortex flow. 

Computed streamlines very near the suction surface ("unwrapped" to lie in a 
plane) and near the endwall are shown in Fig. 2. These limiting streamlines provide 
a clear indication of the strong secondary flow near the endwall toward the suction 
surface and of the formation of the corner vortex (cf., [25], [26]). Additional 
computed results are included here in the form of contour plots of the u, v, 
and w components of velocity (relative to the potential flow direction) and of 
streamwise vorticity in transverse planes at selected streamwise locations. The 
location of the transverse planes, identified by axial step number, is shown in 
Fig. 3. Maximum and minimum contour values are given in Table I. Velocity 
contours at station 2 are shown in Fig. 4. This streamwise location is just 
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downstream of the initial plane, and thus little flow development has taken place. 
The secondary flow proceeds into the endwall boundary layer at the pressure 
surface, toward the suction surface, and out of the endwall boundary layer, but 
the velocities are of small magnitude. At station 12 in Fig. 5, strong crossflow 
in the endwall region and up the suction surface is quite apparent, and crossflow 
velocities on the order of 40% of the maximum streamwise velocity are present. 

The streamwise velocity contours have become distorted due to the corner vortex. 

At station 24 in Fig. 6, the corner vortex fills most of the lower part of the 
passage, and there is considerable distortion of the primary flow. A sequence 
of contour plots of streamwise vorticity £ is shown in Fig. 7. These contours 
also clearly show the formation and development of the corner vortex. 

Finally, although computed results for the static and total pressure are 
ultimately of considerable interest, these are not included here since the 
results were found to contain anomalies which reflect the incomplete state of 
development of the computational algorithm and computer program at the time the 
test case was run. Subsequent detailed examination of the results has led to 
the belief that a straightforward refinement of the computational procedure 
will lead to considerable improvement in the predictions for pressure. Although 
not major, due to the time scales involved it was not possible to develop, 
implement and evaluate the needed changes within the constraints of the present 
study. Despite the inaccurate prediction of pressure, however, the analysis 
in its present state of development provided reasonable results for the secondary 
flow and corner vortex, and appears very promising as a tool for making detailed 
predictions of endwall flow and heat transfer. 
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Fig. 1 


CASCADE GEOMETRY AND 
COMPUTED COORDINATE SYSTEM 
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COMPUTED LIMITING STREAMLINES 


(a) Endwall Surface 



Mid span 





Fig. 3 


LOCATION OF TRANSVERSE PLANES 
AT SELECTED STREAMWISE LOCATIONS 
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TABLE I - MAXIMUM AND MINIMUM CONTOUR VALUES 

(OTHER CONTOURS ARE LINEARLY SPACED) 


Station 

U rain 

u 

max 

v . 
mm 

V 

max 

w . 
mm 

w 

max 

^min 

^max 

2 

0.0 

0.96 

-0.034 

0.014 

-0.008 

0.011 

- 7.7 

- 0.77 

6 

0.0 

1.34 

-0.19 

0.04 

-0.08 

0.15 

- 30.6 

- 2.9 

8 

0.0 

1.77 

-0.41 

0.08 

-0.11 

0.48 

- 59.4 

n 

12 

0.0 

2.25 

-0.87 

0.11 

-0.21 

0.66 

-129.0 

-12.7 

18 

0.0 

2.81 

-0.76 

0.13 

-0.29 

0.86 

-186.0 

-18.4 

20 

0.0 

2.79 

-0.58 

0.13 

-0.32 

0.69 

-158.0 

-15.7 

22 

0.0 

2.71 

-0.45 

0.14 

-0.33 

0.61 

-124.0 

-12.4 

24 

' 0.0 

2.62 

-0.36 

0.16 

-0.34 

0.52 

- 90.1 

- 4.65 
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Pressure Surface 



Suction Surface 



Pressure Surface 


Fig, 4c 



Suction Surface 



Pressure Surface 



Suction Surface 




Fig- 6 


CONTOURS OF VELOCITY 
Station 24 

u - Component v - Component w - Component 
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Pressure Surface 


Fig. 7a 


CONTOURS OF STREAMWISE VORTICITY, C 
Station 2 



Endwall 
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Suction Surface 


Pressure Surface 


Fig. 7b 



Suction Surface 




Fig. 7 


CONTOURS OF STREAMWISE VORTICITY, g 


Station 12 


Station 18 


Station 20 


Pressure Surface 


Fig. 7e 


CONTOURS OF STREAMWISE VORTICITY, £ 


Station 22 


Station 24 


ORIGINAL Pagp ra 

0F p ooR a utU 




Endwall 


10-37 



11.0 DETAILED PROGRAM DESCRIPTIONS 


11.1 MAIN PROGRAM GCCOOR 


Obj ect Compute Coordinates of Circular Arc Duct 

Options IDBG6 = 0 Do not print coordinates 

= 1 Print coordinates 

NDRUM 2 Store coordinates on UNIT NDRUM 


Input Symbols 


DUCTI(l) = di 
DUCTI(2) = d 2 
DUCTI (3) = d 3 
DUCTI (4) = d 4 
DUCTI (5) = d 5 
DUCTI (6) = d 6 
DUCTI (7) = d 7 
DUCTI (8) = d g 


, Duct axial length (ft) (calculated) 
, Duct OD radius (ft) 

, Duct ID radius (ft) 

, Axial location of origin (ft) 

, Radial location of origin (ft) 

, Inlet angle (degrees) 

, Exit angle (degrees) 

, not used 


TITLE (I) 

JL 

KL 

KDS = 1 
IDBGb 
NDRUM 
DDS 


, Title 72 holarith characters 
, Number of streamwise stations 
, Number of streamlines 

9 

, Print option 
, Drum option 

, Mesh distortion parameter 


Output Symbols 


Dl = e 


USC ALE = 

Vr 

R(1,1,J) 

= y id 

R(1,4,J) 

o 

H 

X 

11 

R(2,l, J) 

- y od 

R(2,4, J) 

Q 

o 

X 

II 

JSTEP 


, Roberts Transformation Parameter 
, Metric scale factor 

j> Coordinates of ID radius 

| Coordinates of OD radius 
, Streamwise station number 
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Theory 


The coordinates for a circular arc duct can he obtained from the solution for 
a vortex flow. If W is the complex potential and 2 the complex coordinate, we have 


W = -i/n(z-z 0 ) = s + in 


( 1 ) 


where Z Q is the origin of the circular arc duct and S and n the potential flow 
coordinates. Let 


R = [x-x D ) 2 + (Y-Y 0 ) 2 ] 1/2 


( 2 ) 


6 = tan _1 [(Y-Y o )/(X-X 0 )] 


(3) 


Hence from (1) we have 


S = 9 (4) 

n - -InR (5) 


The potential flow velocity (metric scale coefficient) is given by 

1 

R (6) 

The coordinate functions are then scaled such that 


1 dW 

n v dz 


Hence 


0 < n < 1 

17) 

n = lnCR/R^p) /ln(Rjp/RQp) 

(8) 

s = (Rod “ Rid)*( 0-0 o ) 

(9) 

v = - — - 

(10) 

ln(Rxi>/RoD) R 
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11.2 MAIN PROGRAM GDCOOR 


Obj ect Compute coordinates of Arbitrary Duct 

Options IDBG6 = 0 Do not print coordinates 

= 1 Print coordinates 

NDRUM , Store coordinates on unit NDRUM 


Input Symbols 


JL 


No. streamwise stations 

KL 


No. streamlines 

KDS 

= 1 

No. steps/station 

JLPTS 


No. input data points 

IDBG6 


Print option 

NDRUM 


Drum option 

LFILE 


= 0 do not use 

DDS 


Mesh distortion parameter 

DUCTI(l) 

d l 

Duct length (ft) 

DUCTK.2) 

d 2 

No. knots 

RDII(J) 


Input OD radius (ft) 

RD2I(J) 


Input ID radius (ft) 

Output Symbols 



US CADE 
EPS 


U , Metric normalizing factor 


e , Roberts Transformation Parameter - 
(See Ref. 12 for complete list) 


Theory 

The input data points are smoothed by using a least squares spline fit of the 
second derivative with d 2 knots. This smoothed second derivative is then inte- 
grated to get the smooth contour. The duct coordinates are obtained using the con- 
formal mapping procedure described in Ref. 18. 
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